Relativistic spectral function of nucleon in hot nuclear matter 
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We present a simple calculation of the nucleon self-energy in nuclear matter at finite temperature 
in a relativistic framework, using the real time thermal field theory. The imaginary parts of one- 
loop graphs are identified with discontinuities across the unitary and the Landau cuts. We find that 
in general both the cuts contribute significantly to the spectral function in the region of (virtual) 
' nucleon mass usually considered, even though the unitary cut is ignored in the literature. Also our 

' relativistic spectral function differs from the one in non-relativistic approximation, used in some 

, earlier calculations. 

(N ; 

C^H' Heavy ion collisions provide an opportunity to investigate particle propagation through strongly interacting media. 
^ [ However, only the vector mesons, particularly the p, can at present be studied directly by detecting dileptons, into 
which they decay in the hot, dense media. The media created by these collisions consist, in general, not only of 
00 ' mesons, but also of nucleons. Thus the effects of both mesons and nucleons on the vector meson spectral functions 
C^l . have been extensively studied in the literature [ij. For a more complete picture, the sclf-encrgy of nucleon itself need 

' be investigated 0, Q • The nucleon self-energy function also determines the equation of state of nuclear matter 3 . 
(-H , In this work, we find the one loop corrections to the nucleon propagator at finite temperature and nucleon chemical 
O i' potential in the real time formulation of the thermal field theory [5|. However, to keep an eventual contact with the 
' ] fundamental QCD theory, wc do not start directly with the propagator, namely the two point function of the nucleon 
field. Instead, we consider the same of the nucleon current ri{x) which in vacuum is 

^ '■ r 

' I d^xe'P--'{0\Tr]{x)rj{0)\0). (1) 

<N ■ . . 

^ ' Here ri{x) is built out of three quark fields, so as to have the quantum numbers of the nucleon. We denote its matrix 
CNj ' element between vacuum and nucleon state as 



o 
o 



X 



{0\v{x)\N{p)) ^ Xu{p)e^P-\ (2) 



where u{p) is the Dirac spinor of the nucleon and the parameter A denotes the coupling of ri{x) with the nucleon [8|. 
] In nuclear matter, the vacuum expectation value in Eq. (1) must be replaced by the ensemble average. Denoting 
the time-ordered product of the operators by 0{x), it means 



(0|O(a;)|0) ^ Tr[e-^'^"-''^'^0{x)]/Z, Z = Tr [e-'^^^"^^)] , (3) 

where H and N are the Hamiltonian and the nucleon number operator of the system , /3 is inverse temperature and 
fi is the chemical potential corresponding to N. In the real time version, every two point function, including the 
self-energy we calculate below, assumes the form of a 2 x 2 matrix. But each of these matrices may be diagonalised, 
when it is given essentially by single analytic function, that determines completely the dynamics of the corresponding 
two-point function Q. As this function is simply related to any one, say the 11-componcnt of the matrix, we need 
calculate only this component of the self-energy matrix. 

The 11-component of a free, thermal matrix propagator for a particle is a sum of its vacuum propagator and a 
term depending on the (on-shell) distribution functions of like-particles in the medium through which it propagates. 
The latter term has a universal form, depending only on the bosonic or the fermionic character of the particle [lo| . 
Anticipating the propagators for pion, nucleon and A(1237) to appear in our calculation, we begin by writing their 
11-components. These particles are represented respectively by scalar ((j){x)), Dirac {'ip{x)) and Rarita-Schwinger 
(A^(.t)) fields [llj. For the bosonic propagator, it is 

D^^{ka,k,m^) = A{k,m^) + 2mn{ijj)5{k^ - ml) , (4) 
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FIG. 1: One-loop graphs for the two-point function contributing to self-energy of the nucleon 
where A and n are the vacuum propagator of the scalar field and its equilibrium particle distribution function, 



A(fc,m^) = -2 2 I ■ ' '^H = ^I7— T' w = -y/fc2 + m2 . (5) 

For the fermionic propagators, we introduce 

{Pq,P, to) = A(p, m) - 27riA^(po)'5(p^ - "^^) , (6) 
where N{po) consists of distribution functions for the particle and anti-particle, 

1 



iV(po)=n+(^')^(Po) + "-(^')^(-Po), "±(^0- ^^(^,^^)^^ , = v/p2 + to2 (7) 
in terms of which the 11-components of spin i and | propagators may be written respectively as 



5"(po,p) = i^ + mN)E'\po,p,mN), (8) 

The spectral function of nucleon is obtained from the Dyson equation, giving the complete propagator S' in terms 
of the free propagator S and self-energy S , 

S' ^S- SJ:S' (10) 

where each element is a 2 x 2 matrix in the thermal indices (besides being 4x4 matrices in Dirac space). As already 
stated, they can all be diagonalised to get the respective analytic functions, denoted by a bar, satisfying 

= S -S^S^ (11) 

which may be readily solved as usual. The self-energy function is related to the 11-component of the corresponding 
matrix by , 

ReS(p) = ReS"(p) 

ImS(p) = coth[/3(po -Ai)/2]ImI]"(p) (12) 
The free propagator S turns out to be the same as in vacuum, 

^-^^ p^ -m% + n] ^ ^ 

The calculation simplifies if we take p = . Also restricting to the anti-nucleon pole in Eq.(13), it becomes 

S{po) = ^ -— , (14) 



Decomposing S and S in Dirac space. 



T,s + 7oS.u , s' = S'g + joS'y , (15) 
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FIG. 2: Graph of phase shift for P33 partial wave in nN scattering from our model compared with experiment [2^ (solid circles) 



it follows from Dyson equation that S'^ = S'^. Then letting E = Eg + E„, we get the complete propagator as 



S {Po) = ^ ^ 16) 



giving the spectral function 



— ImE 

^^(^•'^ = j^;;^W^B^WT(h^ ^^^^ 

The graphs which we wish to evaluate are shown in Fig. 1. Because we are interested in finding the nucleon self- 
energy, we retain only the graphs which couple the nucleon current to the nucleon. Also we include A(1237) resonance 
besides the nucleon in the intermediate state. We shall comment later on the contribution of higher resonances. 

The vertices appearing in Fig.l may be obtained from chiral perturbation theory p"3 - [l3 |. The appropriate field 
variable for pion in the effective theory is not 4'{x) introduced earlier, but the SU{2) valued matrix field u{x) related 
to (l)(x) by u{x) = exp(ir • 0/2F^), where Fjr = 93 MeV, the so-called pion decay constant. The effective Lagrangians 
are |15l |. 

C^NN = igAV^^^75^--|^^7M75T>a^0a-f--- (18) 
£^^A = ^iia{u^)l^t''e,a + h.c. = ?^p^^^^-^l+ + ■■■ (19) 

where = i{u^d^u — ud^u^), a,b - ■ ■ = 1,2 and £12 = —£21 = 1, etc. The coupling constants gA and are to be 
determined phcnomcnologically. As is well-known [TtI. [isj. such a model requires form factors at the vertices, which 
we take in the Lorentz invariant form as 

where p and k are the four-momenta of nucleon and pion at the vertices and A is essentially a cut-off on these momenta. 

We first check this model with experimental data on TriV scattering. The interaction (19) allows us to calculate the 
decay width of A + tt as a function of its energy as 



= m^[t ^-(i^)ip r-— - (21) 



Here p is the three-momentum in the centre-of-mass of ttA'' system, 



P = lE^ • ^22) 
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FIG. 3: Branch cuts of self-energy function in po plane for p = 0. 



In this kinematic configuration, the form factor becomes 



The pion-nucleon partial wave / in the P33 channel may now be written in the form 

■^^^^ ^ E^-ml+ imAT{E) ^^^^ 

We take the resonance parameters at the pole position, itia = 1210 MeV and r(mA) — 100 MeV [l^. Taking qa = 2.2 
and A = 400 MeV [l^, we can satisfy Eq.(21) and also achieve reasonable agreement of the phase shift <533 computed 
from Eq.(24) with experiment [2^ (Fig. 2). Also we take qa = 1.26 [2l| and the same form factor at the ttNN vertex. 

Wc now evaluate the self-energies from graphs of Fig. 1. The 11-component of each of the loops has the general 
form 

S"(po,p) = * J ■0^R{p,k)D'\k)E'\p-k) (25) 

where R(j), k) includes the form factors and the factor from interaction Lagrangian at the vertices, as well as the 
spin-dependent factor in the propagator, all of which wc shall write below explicitly. Inserting propagators from 
Eqs.(4) and (6) in Eq.(25), we get 



S"(po,p) 



d^k R{ko) 



(27r)4 (P -ml + ie){(p - fc)2 - m'^ + ie} 
d^k f Nipo-ko)S{ip-k)^-m^) _ n(fco)J((P - m^) 

(27r)3 ^ { fc2 ^ml+ie (p ~ k)^ - + ie 

+^ J ^R{ko)n{ko)N{po - ko)6ik^ - ml)S{{p - kf - m") (26) 

where m denotes mass of baryon (toat or ttia) in the loop and we show the dependence of R on fco only, suppressing 
other variables for brevity. Here the first term refers to vacuum. The second and the third terms arc medium 
dependent, with distribution functions appearing respectively linearly and quadratically. Observe that the third term 
is purely imaginary. 

Carrying out the k^ integration in all the three terms of Eq. (26), we can easily find the real and the imaginary 
parts of E^^. Though ImS^^ contains quadratic terms in distribution functions, ImE, as defined by Eq. (12), turns 
out to be linear in them , 



imi;(po,p) = ^'^y ■ 



d^k 



(27r)34ww' 

[i?(fco = '^){(1 + n{uj) — n^{uj'))S(j)o — lu — uj') — {n{uj) + n_(cj'))(5(po — 1^ + 1^')} 

+i?(fco = -uj){{n{uj) + n+{uj'))5{po + UJ - uj') - (1 + n{uj) - n_{uj')S{po + uj + uj')}] (27) 



where w and lu' are pion and baryon energies. 



U = Jml + fc^ Uj' = ^777,2 + (p - fc)2 (28) 
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FIG. 4: Imaginary (left panel) and real (right panel) parts of self-energy from ttN and ttA loops. Solid curves represent the 
results of our calculation (relativistic, including unitary cuts). Dotted curves result from non-relativistic approximation. 
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FIG. 5: Real part of nucleon self-energy with (solid) and without (dashed) contribution from the unitary cut. 



For the real part of S, we include only the medium dependent pieces, which are given by the second term in Eq. (26) 
alone, 



ReE(po,p) 



R2n + R^n^ Rin — R-sn^ 



(2tt)'^Alull:' \ Po + uj — uj' pq — lu — uj' 
Rin + RjH^ R2n ~ R^n. 
Pq — cj + Lu' Po +UJ + U}' 



(29) 



where Ri, (i = 1, • • • ,4) denote the function R{ko) evaluated at ko = uj, —uj, po — oj' and po + w' respectively. 

In the following we restrict our evaluation to p = 0, when there will be no angular dependence. First consider the 
imaginary part, giving the cut structure [2^[2J]. The delta- functions in the different terms of Eq. (27) control the 
regions of non-vanishing imaginary parts of S, which define the position of the branch cuts. As shown in Fig. 3, the 
first and the fourth terms give rise to the regions, po > {111 + 171^^) and po < —{171 + 111^^) respectively, giving the unitary 
cuts. Similarly the second and the third terms lead to the regions, < Po 5: (™ ~ itt-tt) and — (m — mj^) < Po ^ 0, 
giving the Landau cuts. 

With p — 0, the value of |fe|, fixed by the (5-functions in Eq. (27), is the magnitude of three- momentum in the 
center-of-mass of the pion-baryon system. 



1^,2 ^ {Po - {m + m^)^}{pl - (m - m 



(30) 
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FIG. 6: Nucleon spectral function at different temperatures for a fixed chemical potential (left panel) and at different chemical 
potentials for a fixed temperature (right panel) 



In this frame, the form faetor (23) simplifies to 



Let us define the variables, 



Fm = ^ (31) 



_ pI + ml - _, _ pl-ml+ 

io — , uj — {6Z) 

2pQ 2po 

which actually coincide respectively with w and u' defined by Eq.(28) on all the cuts, except for w on the Landau 
cut, where w = Then the imaginary parts of S for po > are given by the first and the third terms of Eq.(27) as 

ImSfn ) = / ^(^0 = + - n+ioJ')} , on unitary cut .gg. 

Sttpo I ^(^0 = w){n(|w|) + n+(w')} , on Landau cut ^ ' 

The expression (29) for the real part also simplifies for p = to 



ReE(po) = - .^ 2 P / ^4r^ — =2^ — (34) 

where 

h{uj) = uj'{{uj + Ld)Ri - - uj)R2]n - uj{{uj' + uj')Rzn+ - (cj' - aj')i?4?i-} (35) 

Having carried out the evaluation in terms of the (Dirac) matrix- function R{p,k), it remains to write its explicit 
expressions for the two loops, 

Rip,k) = ^(^) F\p,k){2k-plit-k\^ + ^ + mN)}, (ttA^Ioop) (36) 



4 

R{p,k) = U^] F\p,k)(-k'' + ^^^^^^^\i^-l;t + mA), (ttAIoop) (37) 



3 V^^- 

In this calculation, we consider only N and A in the intermediate state, as we expect the contributions of higher 
mass resonances to be small. The reason is that on both the cuts all the distribution functions decrease rapidly 
with the rise of resonance mass, the only exception being n{uj) for pions on the unitary cut, which does the reverse. 
However, unlike the Landau cut, where all higher mass resonances can contribute, the unitary cut, for a fixed upper 
value of the (virtual) nucleon mass, gets contribution from only a finite number of resonances. Thus in the present 
case, where we restrict po ^ 1-5 GeV, only N and A can contribute to the unitary cut. 
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FIG. 7: Comparison of relativistic spectral function (continuous curve) with its non-relativistic limit (dotted curve). 



A question arises in such calculations, whether a non-relativistic approximation could reproduce the relativistic 
results in a quantitative way [2^. To define this approximation, we rewrite E^^{po,p), the spin-independent factor 
in the 11-component of baryon propagator, given by Eq.(6), as 

nil/ „\ ~1 f I — n+{uj') ^ n+{uj') ( 1 — n_(aj') n_(cj') 



luj' po — W -I- ze po — W — le \po -I- w' — ze po + W -I- ze / J 

Here the first two terms describe the propagation of baryon and the last two that of antibaryon. The non-relativistic 
approximation to this propagator consists in retaining only the first two terms above (26l [27| . This approximation in 
turn gives only the first two terms in Eq.(29) for ReS](poiP)- Further we set uj' = m everywhere for baryon [131 • Note 
that we approximate neither the propagator D^^(k) for pion nor its energy-momentum relation. 

In presenting the results of numerical evaluation, we consider, along with our relativistic framework, also the non- 
relativistic one as defined above and another where the unitary cuts are switched off (by omitting the second and 
fourth terms in Eq. (29)), but relativistic otherwise. We note here that the values of and A as well as the fit to 
the experimental data for 633 that we found earlier remain unaffected by the non-relativistic approximation. Fig. 4 
compares the typical behaviour of relativistic results with the non-relativistic ones for the imaginary and real parts 
of self-energy, separately for the two loops - we see that only the real part for ttA loop differ significantly between 
the two. Fig. 5 does this comparison for the real parts between the complete result and the one without the unitary 
cuts, showing significant difference only at higher masses. (The imaginary parts here are, of course, the same as in 
Fig. 4 without the steep lines representing the unitary cuts.) 

Having made these comparisons, we come back to our model in Fig^ 6 to draw the nucleon spectral function at 

different values of T and /i, which are realised in heavy- ion collisions |29l. [30j . As expected, the height of the peak 
decreases with rise of temperature, while it remains about the same within the interval of chemical potential considered 
here [3]| . Finally we again compare in Fig. 7 a typical spectral function of our calculation with its non-relativistic 
limit. 

We also compare our results with two earlier calculations. Leutwyler and Smilga Q use virial expansion to leading 
order to obtain the self-energy for on-shell nucleon in pionic medium. In Fig. 8 we compare their results for the 
imaginary and real parts with those from our model, setting nucleon and A distribution functions to zero. The good 
agreement shows that our model is realistic, if we recall that they evaluate the virial formula with experimental data on 
ttN scattering. Hees and Rapp Q calculate the self-energy in the imaginary time formulation and takes into account 
only the Landau cut for the imaginary part and correspondingly only the first term in Eq. (29) for the real part. 
However, they take higher mass resonances in the loop, validating their model at higher temperatures and chemical 
potentials. 

To conclude, we calculate the self-energy of the nucleon and its spectral function in the real time version of the 
thermal field theory in the relativistic framework. The imaginary part of the self-energy is built out of contributions 
from both Landau and unitary cuts from one loop graphs with ttN and ttA intermediate states. In contrast to results 
in the literature, we find the unitary cut from the ttA loop to contribute significantly in the upper region of (virtual) 
mass of nucleon considered. The nucleon spectral function turns out to be sensitive to non-relativistic approximation, 
establishing the necessity of relativistic treatment for its quantitative determination. 
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